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A MICROCOMPUTER-BASED SAMPLING DIGITAL VOLTMETER 


INTRODUCTION 


Acoustical measurement systems at the Underwater Sound Reference 
Detachment (USRD) of the Naval Research Laboratory need a specialized 
multiple-channel voltmeter with the capability to measure voltage, relative 
phase between channels, and harmonic distortion for pulsed sinusoids at 
frequencies between 0.1 Hz and 100.0 kHz. To answer this need, a sampling 
digital voltmeter was developed around an LSI~-11/23 microcomputer, using three 
analog-to-digital (A/D) converters and an IEEE-488 interface to serve as the 
coumunications link with any appropriate host computer. Using this voltmeter, 
a “simple” computer-controlled measurement system can then be configured with 
the following auxiliary equipment as shown in Fig. 1: synthesizers to 
generate a sinusoidal test signal and a coherent sampling signal, and 
programmable amplifiers or attenuators as needed. (The voltmeter has no 
intrinsic autoranging capability.) 


HOST 
COMPUTER 


ie 
wv 


\EEE-488 BUS 


EXTERNAL ~ SAMPLING — 
are SUL SAMPLE RATE DIGITAL 
GENERATOR VOLUARETER 


(SEE FIG. 2) 


REFERENCE 
FREQUENCY SIGNAL 


POWER 
SYSTEM 
UNOER TEST 


AMP 
: Fig. 1 - Typical measurement system using the 
microcomputer-based sampling voltmeter. 
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APPLICATIONS 


The two immediate applications for such a voltmeter were for the USRD's 
Anechoic Tank Facility (ATF) and Low-Frequency Facility (LFF) measurement 
systems. 


The ATF high-power measurement system requires a three-channel voltmeter 
to measure simultaneous projector drive voltage, drive current, and receiver 
output (hydrophone voltage). Measurements of relative phase between voltage 
and current and the harmonic distortion in all three channels are necessary. 
The practical frequency range is 2 to 100 kHz; the mode of operation is 
usually pulsed. 


The LFF system requires a two-channel voltmeter to measure the 
simultaneous voltage outputs of two hydrophones, generally a probe standard 
and an unknown hydrophone. Relative acoustic phase is required, and harmonic 
distortion is a useful parameter to monitor. The practical frequency range is 
0.1 Hz to 4 kHz; the mode of operation is continuous wave (cw). 


The voltmeter was, therefore, conceived to be a device independent of 
(and transportable between) specific measurement systems at the USRD, with the 
possible exception being that the A/D converters might vary with application 
(e.g., a need to operate at higher sampling rates) and are otherwise 
functionally identical plug-in modules. To get the widest possible range of 
potential host computers, communication through the I[EEE-488 bus was chosen. 
Another requirement was that all external hardware, including the sample-rate 
generator, be controlled by the host computer for maximum flexibility. 
Therefore, all external functions, including automatic ranging of the 
voltmeter, are controlled by the host. 


Since initiating and implementing this signal processing sampling 
voltmeter, two other similar instruments have been described or released: the 
low-frequency sampling voltmeter by the National Bureau of Standards ({1] and 
the sampling network analyzer by the Dranetz Company (2]. The three 
instruments share the sampling ideas but differ in the features and algorithms 
implemented. 


THEORY OF OPERATION 


The voltmeter consists of a microcomputer containing up to three A/D 
converters with direct memory access (DMA) capability (Fig. 2). This 
configuracion allows the computer to digitize up to three input waveforms 
simultaneously, sampling them when triggered by a common external sample-rate 
generator. Where simultaneous measurement of all channels is not desired, 
separate or delayed sampling signals can be given to individual A/D 
converters. The present voltmeter is designed to digitize three waveforms of 
up to 1024 samples each; these buffers could be expanded to 3900 samples (see 
Appendix F). 


Fig. 2 - Block diagram of the microprocessor- 
based sampling voltmeter. 


The DMA interface loads 12-bit digitized values directly into memory at 
rates up to 100 kHz per channel. The individual sampled waveforms can be 
transmitted to the host for processing; alternately, the slave can compute 
sinusoidal and statistical parameters from these waveforms and return them to 
the host. Parameters are computed for the fundamental frequency. (Throughout 
this report “fundamental” frequency refers to the sinusoidal test frequency 
being be measured, no matter how many cycles of this frequency are included in 
a measured sequence. ) 


The sinusoidal quantities amplitude, phase, and harmonic distortion are 
computed for each channel using single frequency discrete Fourier transforms 
(SFDFT). The statistical quantities average, standard deviation, and the 
positive and negative extrema are obtained for each channel to characterize 
noisy or nonsinusoidal signals. In voltmeter terms, these statistical 
parameters are the dc, true ac rms, and positive and negative peak voltages. 
To measure the coherence or correlation between channels, the second order 
joint moment statistic, E(xy), can be computed. For the case of the two 
signals, x and y; proportional to current and voltage, the joint moment is 
proportional to power. The covariance and cross correlation can be readily 
obcained from these parameters. 


COMPUTATION PROCEDURES 
Sampling Criteria 


In the usual case of interest (pulsed sinusoids), sampling errors are 
minimized through coherent sampling, which requires an integer number of both 
test frequency cycles and sampling periods: 


where fp is test frequency 
Es is sampling frequency 
n is samples per sequence (integer) 
m {8s cycles per sequence (integer). 


Lf the sampling frequency, 


f, 75 f , (2) 


is greater than twice the test frequency, the Nyquist theorem is satisfied, 
and satisfactory results will be obtained. If the necessary sampling rate is 
too fast for the A/D converter to operate properly, undersampling can be used; 
by choosing fewer than two samples per cycle while retaining an integer number 
of samples and cycles per sequence, an effective higher sampling rate is 


attained. A comparison of the oversampled waveform shown in Fig. 3 with the 
undersampled waveform in Fig. 4 clearly demonstrates how undersampling can 
synthesize a waveform (Fig. 5) effectively sampled at a higher rate, although 
with fewer cycles.overall. Undersampling will only work properly on periodic 
waveforms, however. Computed data in Table 1 demonstrates that peak voltage 
and harmonic distortion may be affected slightly by undersampling. 


iN 


Fig. 3 - Steady-state sinusoidal waveform sampled for 16 cycles at 17 
samples per cycle: 272 samples (oversampling). 


AN 


Fig. 4 - Steady-state sinusoidal waveform sampled for 16 cycles at 1.0625 
samples per cycle: 17 samples (undersampling). 


Fig. 5 - Steady-state sinusoidal waveform synthesized from the 17 samples 
shown in Fig. 5. 


Table 1 - Comparison of oversampled and 
undersampled measurements. 


OVER UNDER 
COMPUTED PARAMETER SAMPLING SAMPLING 

RMS VOLTAGE (STD DEV) 1.047 

PEAK VOLTAGE .485 

DC VOLTAGE (MEAN) .013 

TOTAL POWER . 102 

RMS’ VOLTAGE (DFT) 

RELATIVE PHASE 

HARMONIC DISTORTION 


ST II ETE PTET SOE 0 TET OC SET TAREE 
Normal: 2 cycles ac 26 samples per cycle (52 samples 
total). 
Unders@mpled: 25 cycles st 1.04 samples per cycle (52 
samples total). 


Sinusoidal Parameters 


For sinusoidal testing, the signals are of the form 


y(t) = 2a cos(2nfit +o) 


g 
+ Y2A,_cos(2mkf,_t +o, ) #n(t) , 
a An T km 


where the terms are, respectively, the fundamental, the harmonics, and the 
additive noise. 


Least squares estimates of the rms amplitude, A,, and phase, a? of the 
fundamencal are obtained from 


P 1 
7 2 2 
A, = 7% [Re2(x,) + Im*(y,)] f 


7 ae [ macx,>/Recy,)| ; 


l n-1 
oe = 2, y, exp(~j2rim/n) ’ 


which is the mth component of the DFT of the n-point sequence y, containing m 
cycles of the fundamental; that is, the single frequency discrete Fourier 
transform (SFDFT). With coherent sampling and in the absence of noise, An 
and >, are exact. 


Harmonic distortion, D, is defined in terms of the ratio of power in the 
harmonics to power in the fundamental. 


oof hating?) Me. 


k=2 


Since the sequences are relatively short and few lines are required, 
Goertzel's algorithm [3,4] is an efficient implementation of the SFDFT. 


Digital Filtering 


The evaluation of the DFT at a single frequenc, as in Eq. (6) can be 
interpreted as a digital filter (5] with a corresponding enhanced signal-to- 
noise ratio. However, the filter quality, the skirts and side lobes, can be 
significantly improved by weighting or windowing the data sequence and then 
doing the SFDFT. A particularly effective choice of data weighting sequence 
is the Kaiser I,-sinh window [6]. An example of a weighted data sequence is 


shown in Fig. 6, and the corresponding digital filter responses are shown in 
Fig. 7. 


Fig. 6 - Sinusoidal waveform of 10 cycles (product 
with Kaiser window with 40 dB of spectral side 
lobe attenuation). 
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Fig. 7 ~ Comparison of the filter effect of an unweighted SFDFT 
(solid line) with thace of a SFDFT with the Kaiser window with 
40 dB of side lobe attenuation (dashed line). The frequency 
axis is normalized to the center frequency f, which 
corresponds to the mth spectral line. The data sequence 
includes 10 cycles. 


The filters are centered on the mth spectral line or test frequency, fr: 
Following convention, the filter bandwidth is defined as the main lobe width 
or the distance between the minima adjacent to the mth line. For the 
unweighted SFDFT, the fractional bandwidth or main lobe width is 


oe “ 
en Rectangular 


where m is the number of cycles of fy, in the data sequence. 


, (8) 


B\(N 


The reduced side lobes of the Kaiser filter are achieved at the expense 
of an increased main lobe width 


Af 2} 6(R + 12) 
jel e({l} 2’ , (9) 
(ae (3) 1? 


where R is the attenuation in dB of the first side lobe. Figure 8 illustrates 
the increase in filter main lobe width as the side lobes are reduced. 


ATTENUATION {dB} 


I/tm (FREQUENCY NORMALIZED 10 THE CENTER FREQUENCY Im) 


Fig. 8 - Comparison of filter main lobes for unweighted 
SFDFT and Kaiser-weighted SFDFT with 40 to 100 dB 
side lobes. The side lobes are not plotted for 
clarity. The data sequence includes 10 cycles (m =* 10). 


The data window or weighting sequence, W;, 1s computed in the host and 
transferred to the voltmeter where the window is applied in estimating 
amplitude and phase using the SFDFT of the weighted sequence 


a-l 
De wy y,exp(-j2nim/n) . 
i290 


The amplitude estimate, A,, is corrected by the window scale factor 


Table 2 illustrates that the Kaiser window does not alter or bias the 
amplitude and phase estimates at the test frequency. However, the harmonic 
distortion, particularly for low distortion, is biased by the window; larger 


sidelobe attenuations reduce this bias by reducing the amplitude of the 
harmonics created by the window. 


Table 2 - Amplitude, phase and harmonic distortion 
estimates with and without a Kaiser window. 


eS fg Sr sD ES SS Se ST 


SIMULATED DATA 


MEASURED DATA 


t A 100D PERCENT simcure DISTORTION 
a a PERCENT HABMONIC 
WINDOW RMS AMPLITUDE PHASE DEGREES DISTORTION of 0.1% 1% 

No Window 1.054 -123.3 0.060 0.000 9.100 1.000 
- 30 dB Side Lobe 1.054 -123.3 0.071 0.032 0.073 0.972 
- 40 dB Side Lobe 1.054 -123.3 0.077 0.018 0.084 0.984 
- 50 dB Stde Lobe 1.054 -123.3 0.082 0.009 0,092 0.992 
- 60 dB Side Lobe 1.054 ~123.3 0.087 0.004 0.097 0.997 
- 70 dB Side Lobe 1.054 -123.3 0.097 0.002 0.099 0.999 
- 80 dB Side Lobe 1.054 -123.3 0.093 0.001 0.099 0.999 
- 90 4B Side Lobe 1,054 -123.3 0.095 0.000 0.100 1.000 
-100 dB Side Lobe 1.054 -123.3 0.097 0.000 0.100 1.000 


Computations performed on 20 cycles of data at 32 points per cycle. Measured data was the output of 
a synthesizer: simulated data was a computed sinewave with distortion added to the second barmontc as 
labeled. 


When long sequences or long integration times are used to combat noise, 
the computation time can be significantly reduced by averaging. This 
efficiency is possible if the total sequence consists of two or more 
subsequences with an integer number of samples and cycles per subsequence. 
The simplest example is a date sequence of m cycles of the fundamental test 
sinusoid with each cycle containing n samples--this data sequence can be 
compressed from m cycles to one cycle by averaging. Both the SFDFT and 
windowed SFDFT computations can be speeded up through averaging. (The window 
is applied before averaging.) 


Statistical Parameters 
For signals which may be arbitrary periodic waveforms or noise-like data, 

statistical parameters are more appropriate in characterizing the data than 

are the SFDFT derived amplitude and phase (although naturally these parameters 


can be computed for sinusoidal signals as well). 


The voltage extremes or range of the waveform are obtained by inspection 
of the data sequence, and both the maximum and minimum voltages are found. 


The average or dc voltage is computed with 


a 9. te 
Yao 7% 77 hy Ms * (12) 


(13) 
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The choice of factor 1l/n instead of 1/(n-1) is to satisfy the Parseval 
relationship so that the rms amplitude of an ideal sine wave is the same 
whether estimated with the SFDFT or with the true rms methods. 


The second joint moment ° 


1 2 
E(xy). = = & iY, 9 (14) 


is analogous to power if the data sequences X; and y; are proportional to 
current and voltage. Power measurement was the motivation for including this 


joint moment estimate. The total ac power, P,.., is computed using 


Pac ™ E(xy) - xy - (15) 


All of the summations in Eqs. (10) through (15) can be interpreted as 
Euler-Maclaurin approximations to integrals of the continuous time 
functions. For functions that can be represented as a finite number of 
harmonically-related sinusoids, this simple approximation is exact [7]. 


SOFTWARE FUNCTIONS 


The voltmeter software is organized around a monitor routine awaiting 
input from the host via the [EEE-488 interface (as shown in Fig. 8). The 
input consists of ASCII strings of variable length--the first two characters 
being a code for setup or control. Strings are terminated by an 
end-or-identify (EOL) sent with the last character. 


“Setup” codes are accompanied by an encoded byte string specifying the 
value foc a certain parameter, such as the programmable gain for a specific 
A/D converter. Setup strings tell the program what values are to be used, but 
do not perform any operation such as physically setting the gain of an A/D 
converter. No default parameters are present within the program; all 
parameters must be specifically loaded. 


“Control” codes (two-byte strings) command the voltmeter to perform 
certain functions; i.e., make an A/D conversion, compute certain parameters, 
transfer data to the host, etc. 


The voltmeter can be requested to compute and return various parameters 
as requested, or to return the original sampled waveform to the host for 
special processing not supported by the voltmeter (plotting, for example). It 
is necessary to point out that although commands are transferred to the 
voltmeter as ASCII byte strings, data is received from the voltmeter in PDP-11 
format (two bytes forming the l6-bit integers for sampled data, or four bytes 
floating point variables for all computed parameters). This was done in the 
interest of faster data transfer and maintaining original data precision for 
NRL-USRD applications. This feature prevents full compatibility with non-PDP 
11 host processors, although subroutines to encode integer and floating-point 
variables exist as patt of the display functions, and the program could be 
modified accordingly. 
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Fig. 9 - Flow chart of the basic monitor 
logic in the voltmeter software. 


A special control code (VM) allows the voltmeter to measure, compute, and 
display as its own pace without host control; the host regains control by 
generating a group-execute-trigger (GET) on the IEEE-488 bus. The voltmeter 
displays data on either the console terminal or a 24-character plasma 
display. Display styles allowed are: count of A/D conversions and errors; 
voltage, phase, and distortion of one specific channel; peak voltage of all 
‘hree channels; or voltages of two channels and the phase between the two. 


The computation commands follow this sequence: the buffer of 12-bit 
digitized integer values is converted to floating-point numbers representing 
voltages within the A/D converter range (+10 V); a window function is applied, 
1£ desired; the data sequence in the buffer is averaged into a smaller 
subsequence, if desired; the specific computation is made; the results are 
scaled to account for the A/D converter gain and corrected for amplitude 
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change due to the window function, if used; and the data are transmitted to 
the host. The use of windowing and subsequence averaging seems not very 
useful for the computation of nonsinusoidal parameters. 


The four primary commands are: “D” (compute voltage, phase, and 
distortion from SFDFT); its subset, "C” (the same without distortion), used 
for more rapid computation where distortion is not needed; "E” (total rms, 
peak, and de voltage); and a similar subset, “F" (absolute peak only), a rapid 
tool for system autoranging or overload detection. Relative computation times 
are shown in Table 3. Note that these times refer to a PDP-11/23 computer 
with both the KEF-11 floating-point processor and the new FPF-11 floating- 
point hardware. The FPF~11l offers as much as a factor of 5 improvement in 
computation speed. 


Table 3 - Sample computations measured for a PDP-11/23 
with the KEF~-11 and FPF-L1 floating-point processors. 


COMPUTATION TIME IN SECONDS 
NO SUBSEQUENCE AVERAGING WITH SUBSEQUENCE AVERAGING 


RO WINDOW WITH WINDOW NO WINDOW WITH WINDOW 


COMPUTATION STYLE KEF-11 Ypr-1i1 KEF-11 FPF-11 KEF-11 FPF-11 KEF-11 FPF-11 


SINUSOIDAL PARAMETERS 
Voltage, Phase 0.494 0.149 0.210 0.098 
Voltage, Phase, Distortion 2.627 0.506 0.383 0.133 


NONSINUSOIDAL PARAMETERS 
rms, Perak, dc Voltage 0.350 0.182 0.198 0.100 
Peak Voltage 0.189 0.096 - - - a 
Total Power 1.488 0.477 - - 0.927 0.373 


ea ST ET 


Average of 100 measurements of 20 cycles of data at 40 points per cycle (800 points). Distortion computed from 
first 7 harmonics. Subsequence averaging to 1 cycle. Times include data transfer on IEEE-488 bus (about 
0.016 second). 


Collateral computation commands include the ability to compute the total 
rms voltage and phase of a specific harmonic and to compute the power 
parameters second joint moment (ac and dc power) and covariance (ac power). 
The power parameters are the total power generated by all harmonics, and it is 
assumed that the host will scale the data properly to account for the 
calibration of the current measuring circuitry. 


Further details of the software are to be found in the appendices and the 
liberally-commented source code located there. Appendix C lists various error 
and status messages; Appendix D gives listings for sample host routines to 
drive the voltmeter; and Appendix E includes complete listings of all software 
modules loaded into the voltmeter, with external discussion where appropriate. 


SUMMARY 


A multiple-channel microprocessor-based sampling voltmeter has been 
developed for use in the USRD measurement systems. This voltmeter is capable 
of measuring amplitude, phase, and distortion, and of using several 
computational algorithms to derive these results as the application 
requires. The original voltmeter has been used successfully in the Anechoic 


Tank Facility as a high-voltage impedance measurement system, and this more 
powerful version has been greatly expanded to be used in Low-Frequency 
Facility, or other measurement systems as needed. 
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APPENDIX A 


Voltmeter Hardware 


The hardware chosen as the heart of the voltmeter was the Digital 
Equipment Corporation (DEC) PDP-11/2 microcomputer (later upgraded to the PDP=- 
11/23 with a KEF-1l floating-point processor), because it was a device that 
the USRD had both the experience and facilties to handle. 


For A/D converters, ADAC Corp. models with DMA interfaces were chosen: 
model 1012, whose maximum sampling rate of 35 kHz was satisfactory for the LFF 


voltmeter, and model 1012AD, whose maximum sampling rate of 100 kHz was 
necessary for the ATF voltmeter. 


Other interface modules used were: 


MXVLI1-A (DEC) Multifunction module: 
2-serial line units 
16K words RAM memory 
4K words ROM memory (Intel 2732 EPROMs) 


GPIBVL1-L (National Instruments) [EEE-488 bus interface 


DRV1LL (DEC) L6—bit parallel 1/0 card 


The parallel 1/0 card is used only in pulsed systems to synchronize the 
measurements with the beginning of a toneburst. 
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APPENDIX B 


Voltmeter Commands 


SETUP COMMANDS (# is channel number): These consist of a 2=byte code and a 


AI 


AN 


AV 


CT 


OD 


DS 


FC 


PA 


PR 


decimal integer string for the parameter value. 


Number of times to automatically retry a measurement if an A/D 
error occurs 


Which A/D-DMA device generates an interrupt when the data transfer 
is complete (1, 2, 3) 


Which A/D converters are to be used; additive, 
where 1 = #1, 2 = #2, 4 = #3, (7 = all three) 


Averaging mode (0 = no averaging, >O = number of cycles in the 
averaged subsequence) 


Computation type (for future use) 

Display device: serial line 1 (console) or 2 (plasma display) 
Display style: 0 (conversion count), 1-3 (voltage, phase, 
distortion of channels 1-3), 4 (peak voltage in all 3 channels), 
5 (voltages in channels 1-2 and phase between) 

Sample (clock) frequency in Hz 


Signal frequency in Hz 


Gated system (0 = not gated, 1 = use parallel 1/0 card to 
synchronize with toneburst) 


Programmable gain for A/D converter (1, 2, 5, 10) 

Number of harmonics used to compute distortion 

Which harmonic to compute (see B# - CONTROL COMMANDS) 

Transfer message to display, a byte string of up to 40 characters 
Number of cycles per sequence 


Number of points to convert per sequence 


Number of points per sequence to return to host 


PS 


P# 


CONTROL 


BH 


CW 
Ci 


Di# 


R# 


VM 


WI 


x2 = Pee a ce NNN tere SAREE 2 easy 


First point to return to host. In practice, the first few 
converted points are bad, therefore, PA = PR + PS points are 
converted and PR points are computed 

Port (channel) for A/D Multiplexor (0-15) 


Window mode (0 = don't use window, 1 = use window) 


COMMANDS (2-byte control codes; # is channel number) 


Compute and return voltage and phase of the harmonic defined by HC 
parameter 


Compute a sine wave for test purposes 
Compute and return voltage and phase of the fundamental 


Compute and return voltage, phase, and distortion of the 
fundamental 


Compute and return statistical data (rms, dc, positive and 
negative voltage extremes) 


Compute and return absolute peak voltage 

List setup parameters on the display 

Perform A/D conversion and return status 

Compute and return total power (second joint moment) 
Transfer setup parameters to host 

Clear conversion and error counters 

Read back a data sequence 

Read back the average of a window and the weighted sequence 


Compute and return total power (second joint moment) and 
covariance (ac power) 


Measure the time used in data transfer over the I[EEE-488 bus 
Set to free-running display voltmeter mode 


Prepare voltmeter to read a computed window sequence from the host 
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SPECIAL LEEE-488 COMMANDS 


. EOL End or identity - must coincide with last character sent by the 
, host 
GET Group execute trigger - command will interrupt whatever the 


voltmeter is doing and reinitialize 


APPENDIX C 
Voltmeter Status Messages 
Status is returned as an integer variable that is the voltmeter's 


acknowledgment of a measure command. It refers to either the legality of the 
command or the success of executing it. 


STATUS MEANING 
1 Normal return 
~l A/D error in channel 1 
~2 A/D error in channel 2 
-4 A/D error in channel 3 
~10 DMA error in channel 1 
~20 DMA error in channel 2 
-40 DMA error in channel 3 
-100 A/D setup error 
-101 Too many points for program buffer 


Errors -1 through -40 are added together in case of wultiple errors, and 
thus cover all numbers between ~1 and -77. 


The voltmeter processor will halt for the following severe errors: 


ADDRESS REASON FOR HALT 
002036 Floating-point exception (can proceed) 
002054 Function (sine, atan, sqrt) error (can proceed) 
003370 Too many samples to compute (can proceed) 
006224 Illegal display mode (can proceed) 
013236 IEEE~488 interface illegally strapped (fatal) 
21 
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APPENDIX D 


VYoltmeter Control Subroutines 


The following describes the protocol necessary and the available general- 
purpose subroutines for controlling this sampling voltmeter from a host 
processor. The National Instruments [EEE-488 bus driver and its software 
interface IBUP are a prerequisite [8]- 


Listings for the following subroutines (written at NRL-USRD) are included 
in this appendix. 


1. Primitive subroutines 


RDLSI Poll the bus and if successful 
read a byte string from the voltmeter 


WRLSL Write a byte string to the voltmeter 


2. Higher-level I/O subroutines 


WRLPA Send an integer parameter to voltmeter 
WRLVM Send an integer array to voltmeter 
WRMEA Send measurement command to voltmeter, 


wait an appropriate moment. 

enable specific clock synthesizer, 
receive status from voltmeter, 
disable specific clock synthesizer 


3. Example of voltmeter 1/0 


METST Test routine to command a measurement, 
read back a data sequence, and 
convert to volts 
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RDLS! . FTN; 161 


000! 


Cc 
Cc 
Cc 
Cc 
Cc 
Cc 
Cc 
Cc 
Cc 
Cc 
Cc 
i Cc 
' 0002 
: 6003 
' Cc 
; 0004 10® 
Cc 
; Cc 
j Cc 
: cD 
cDdi1e 
c 
Cc 
Cc 
0005 
0006 
0007 
D 
D120 


0009 68 
6010 
Oe 
i" D 
D210 
D 


4012 


9013 


0014 


8015 


PROGRAM 


Name 


SCODE1 


FORTRAN [V-PLUS 


SECTIONS 


Size 


600216 
SPDATA 0600010 


3IDATA 000040 
S3VARS 900004 
Total Space Allocated = 900272 


No FPP [nstructions Generated 


-@ 21:07:42 25-Mar-82 


7TR: BLOCKS/WR 
SUBROUTINE RDLSI (IADR, IBUF,ILEN, ISW) 
SUBROUTINE TO READ A SYTE STRING FROM THE LSI VOLTMETER. 


V3 


ARGUMENTS : 
IADR = GPIB BUS ADDRESS 
IBUF = BYTE-BUFFER TO LOAD WITH STRING 
{LEN = LENGTH OF STRING EXPECTED 
iSW = STATUS OF READ PROCESS 


ISW = 1 (SUCCESSFUL) 
AUTHOR: R. E. SCOTT, JR. (NRL’USRD); VERSION: JANUARY 1982 


BYTE IBUF({LEN) 
INTEGER GPIB 


CONT I NUE 


TEST FGR SRQ (CURRENTLY NOT WORKING WITH LSI) 
J=GPIB(13) 

IF (J.NE.1) WRITE (6,110) J 

FORMAT (' RDLSI TEST-SRQ FAILURE: '14) 

IF (J.NE.1) GOTO 100 


POLL FOR RESPONSE 

J=1BUP (6, IADR) 

IF (J.EQ."102) COTO 20@ !SUCCESS 
IF (J.EQ@Q.-6) COTO 236 {TIMEOUT 
WRITE (6,120) J 

FORMAT (’ RDLSI POLL FAILURE: ' 14) 
GOTO 100 {TRY AGAIN 


READ THE BYTE STRING 
J=IBUP (1, 1ADAR, I BUF, ILEN) 

IF (J.EQ.ILEN) ISW=1! ! SUCCESS 

IF (J.NE.ELEN) ISW=J 1 FAILURE 

WRITE (6,210) J,1SW 

FORMAT (7° RDLSI STATUS: J='13° ISW='13) 
WRITE (6,220) IBUF 

FORMAT (" RDLSI STRING=' 1204) 


RETURN 


CONTINUE 
WRITE (6,240) J 

FORMAT (° RDLSI TIMEOOT: J='13). 
RETURN 


END 


Attributes 


RW,1,CON,LCL 


4 RW.D,CON,LCL 
16 RW,D.CON,LCL 
2 RW,D,CON,LEL 
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FORTRAN [{V-PLUS V3.0 21:07:38 25-Mar-82 Page ! 


WRLSI . FTN; 1@1 7TR: BLOCKS/WR 
000! SUBROUTINE WRLSI (IADR,STRING, ICNT,ISW) 
Cc 
Cc SUBROUTINE TO WRITE A SYTE STRING TO THE LSI VOLTMETER. 
Cc 
Cc ARGUMENTS : 
c IADR -- BUS ADDRESS OF DEVICE 
Cc STRING-- BYTE ARRAY TO BE SENT 
Cc {CNT -- LENGTH OF ARRAY TO BE SENT 
c isw -~ STATUS OF OPERATION 
Cc ISWs1 1S SUCCESSFUL 
Cc IGW=* UNSUCCESSFUL (* = GPIB ERROR) 
Cc 
Cc AUTHOR: R. E. SCOTT, JR. (NRLYUSRD); VERS(ON: JANUARY 1982 
Cc J 
6002 BYTE STRINGCICNT) 
Cc 
c WRITE THE STRING 
9003 J=1BUP(@, IADR, STRING, ICNT) 
Cc 
C DID [T WORK? (IF [SW = ICNT, SUCCESS;. 
c iF NOT. ISW {3 GPIB ERROR STATUS. 
0004 IF (J.NE.ICNT) ISW3J 
6065 iF (J.EQ.1CNT) [SW=1 
Cc 
dD WRITE (6,90) STRING 
D9 FORMAT ¢' WRLSI STRING=° <ICNTDAIL) 
5 WRITE (6,100) I[ADR.J,1SW 
D1e6e FORMAT (° WRLSI STATUS: ADDR2'03° J='13° [SWs' 13) 
Cc 
0006 RETURN 
0007 END 


PROGRAM SECTIONS 


Name Size Attributes 
SCODE! 900134 46 RW, 1,CON,LCU 
SPDATA 060004 2 RW,D,CON.LCI. 
SIDATA 0090932 13 RW,D,CON,LCL 
SVARS 000062 ! RW,.D,CON,LCL 


Total Space Allocated = 000174 62 


No FPP Instructions Generated 


FORTRAN IV-PLUS V3.0 21:08:16 25-Mar-82 Page 1 
WRLPA .FTN; 161 


0601 


0002 
9003 


0004 
0005 


0006 
0007 


0008 
9009 
0010 


O@11 
O012 


013 
C014 


PROCRAM 
Name 


SCODEI 
SPDATA 
SIDATA 
SVARS 


AANAAARAMAAAMAA 


QO 


[<) 
S 
.) 


oo 


31® 


AAGAAAAAAAGSA 


7TR: BLOCKS/WR 
SUBROUTINE WRLPA (IADR, ICODE,1VALUE, ISW) 


SUBROUTINE TO WRITE A CODE AND AN INTZCER VALUE TO 
THE LS!I-VOLTMETER. AND READ AN ACKNOWLEDGEMENT . 
1ADR = GPIB BUS ADDRESS OF VOLTMETER 
ICODE = TWO-BYTE PARAMETER SETUP CODE 
IVALUE = POSITIVE INTEGER PARAMETER TO LOAD 
ISw = STATUS OF OPERATION 
ISw = 1 (SUCCESSFUL AS ACKNOWLEDCED) 
{SW => * (FAILURE * AS ACKNOWLEDGED) 


AUTHUR: R. E. SCOTT, JR. (NRL’“USRD); VERSION: JANUARY 1982 


BYTE STRING(8) 
BYTE ICODE(2) 


INCORPORATE THE CODE INTO THE STRING 
STRING(: 3 = (CODEC!) 
STRING(2)= [CODE(2) 


ENCODE PARAMETER INTO SIX-CHARACTER ASCII STRING 
ENCODE 1(6,160,STRING(3)) [VALUE 
FORMAT (16) 


REPLACE LEADING SPACES WITH ZEROS 
DO 110 [23,8 
IF (STRING(1]).£Q.° ') STRING(I)='@° 


NB=8 


WRITE THE STRING TO THE CPIB-BUS 

WRITE (6,200) IVALUE,STRING 

FORMAT (' WRLPA VALUE='16’ STRING="8A1) 

CALL WRLSI (IADR,STRING,NB, [SW) 

IF CISW.NE.1) CALL TSTERR (6, 1SW) TERROR IS ??7 

IF cISW.EQ@.-6) J=IBUP(1@) !1TIMEOUT CLEANUP 
WRITE (6,300) ISW 

FORMAT (' WRLPA (WRITE) STATUS: ISW='13) 


READ REPLY IF GIVEN 

CALL ADLSI (IADR, ISTAT,2,1SW) 

IF (iSW.NE.1) CALL TSTERR '6,1SW) !ERROR [S 727? 

IF ctSW.EQ.-6) J=1BUP(1@) !TIMEOUT CLEANUP 
WRITE (6,510) ISW 

FORMAT (' WRLPA (READ) STATUS=° 13) 


SUBSTITUTE ACKNOWLEDGED STATUS FOR GPIB STATUS 
(SW=iISTAT 

RETURN 

END 


SECTIONS 


Size 


000252 
960610 
000032 
000014 


Attributes 


85 RW,1,CON,LCL 
4 RW.D,CON,LCL 
13 RW ,D,CON,LCL 
6 RW.D, CON, LCL 


Total Space Allocated = 900990 108 


No FPP Inetructions Generated 
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FORTRAN IV-PLUS V3.0 21:08:36 25-Mar-82 Page |! 
WRLVM. FTN; 161 


6001 


6002 
0003 
0004 


0005 
0066 


807 
0008 


600? 


0010 
091i 


eO12 
9013 


PROGRAM 
Name 


SCODE! 
S3PDATA 
SIDATA 
SVARS 

STEMFS 


ANAAARAGAAAAAAG 


Aaa 


loka 
~) 


QOADAAM& 
S 


SECTIONS 
Size 


9OO254 
600624 
0098652 
9000390 
0900062 


7TR: BLOCKS’ WR 
SUBROUTINE WRLVM CIADR, ICODE,1IVAL, ICNT,1SW) 


SUBROUTINE TO WRITE A BYTE STRING TO THE LSi-VOLTMETER 
PREFIXED BY A TWO-BYTE PARAMETER CODE. 


TADR = GPIB BUS ADDRESS 

ICODE = TWO-BYTE PARAMETER CODE 

IVAL = ARRAY TO SEND TO VOLTMETER 

ICNT = NUMBER OF INTEGERS TO SEND 

{sw = STATUS OF ACKNOWLEDGEMENT 
A 


IVAL COULD BE 4-BYTE REAL VARIABLE (ICNT=2). 
IVAL COULD BE NOTHING (ICNT=@; SEND CODE ONLY). 


AUTHOR: R. E. SCOTT, JR. (NRL/US"N); VERSION: JAZUARY 1982 


INTEGER ICODE 
INTEGER IVAL‘8&; 
UNTEGER STRING(1@) 


PREFIX CODE iNTO STRIRC 
STRING(1)=i CODE 
(F (ICNT.E@.3) GOTO 66 


ENTER VALUE INTC STRING, TWO BYTES AT A Ti ME 
DO 59 1-1,iCNT 
STRING( 14+ )=1TVALCT) 


COMPUTE STRING LENGTH 
NB=2+2*ICNT 


WRITE THE STRING 

CALL WRLSI (1ADR,STRING,NB, SW) 

IF CISW.NE.1) CALL TSTERR (6, SSW) YERROR iS 2??? 

IF (ISW.EQ.-6) J=1BUP(10) !TIMEOUT CLEANUP | 
IF (C'SW.NE.1) WRITE (6,100) ISW ; 
FORMAT (° WALVM STATUS: [SW='13) 


HEAD REPLY IF GIVEN 

GALL RDLSI CLADR,ISTAT,2,15W) 

IF Ctsw.NF.1) CALL TSTERR (6,1SW) {ERROR 1S 7??? 

IF (ISW.20.-6) J=fBUPCIO) !ITIMEOUT CLEANUP 
i (ISW.NE.1) WAITE (6,110) ISW 

FORMAT (' KDLSI STATUS: i1S%:* 13) 


EXCHANGE GPIB STATUS FOR ACKNOWLEDGEMENT STATUS 
[SW= {STAT 


AETURN 
SND 
Attritbates 
78 KW, 1 ,.CON,LCL 
= RW,D,CON, LCL 
13 RW,0.CON,LCL 
i2 RW,.D,CON,LCL 


i RW,D.CON,LCL 


Total Space Allocated = 990324 106 


No FPP Instructions Generated 
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FORTRAN IV-PLUS V3.6 21:09:27 23-Mar-82 Page |} 


WRMEA . FTN; f 7TR: BLOCKS/ WR 

6001 SUBROUTINE WRMEA (IADR. ICODE, FREQ, ITIME,ISTAT) 
c 
Cc SUBROUTINE TO WRITE A “ME" (MEASURE) TO THE LSI-VOLTMETER, 
Cc (WAIT), ENABLE SAMPLING SYNTHESIZER, RECEIVE ACKNOWLEDCEMT , 
Cc AND DISABLE SAMPLING SYNTHESIZER. 
Cc *%*x SPECIFICALLY FOR LOW FREQUENCY SYSTEM SYNTHESIZER «xx 
C 
Cc IADR...... 1EEE~488 BUS ADDRESS FOR VOLTMETER. 
Cc ICODE..... MEASUREMENT CODE (ME). 
Cc FRE@...... SAMPLING FREQUENCY IN HZ. 
Cc ITIME..... TIME TO WAIT BETWEEN MEASURE COMMAND AND EABLING 
Cc SAMPLE-RATE CENERATOR (MSEC). 
Cc ISTAT..... MEASUREMENT STATUS. 
C 
Cc AUTHOR: AR. E. SCUTT, JR. (NRL“USAD); VERSION: JANUARY 1962 
Cc 


0002 ISTAT=0 !STATUS: NO MEASUREMENT 
C COMMAND 4 MEASUREMENT 
6003 CALL WRLSI (IADR. 'ME’ ,2,1SW) 
004 IFCISW.NE.1)CALL TSTERR (6,1S5W) !BRUS ERROR? 
D [F (iSW.NE.t) WRITE (6,106) ISW 
D:¢0 FORMAT (' WKMEA WRITE-“ME" STATUS: '13) 
0005 IF (CISW.NE.() RETURN 'GIVE UP IF BUS ERROR 
Cc 
0605 CALL WAIT CITIME,1,K) !WAIT FOR SLAVE TO BE READY 
9007 CALL SM102 (°24,FREQ,.'E °,JSW,ISW)!tENABLE SYNTHESIZER (CLOCK) 
0008 CALL RDLSI CLADR,ISTAT,2,1SW) !GET STATUS FROM SLAVE 
0009 IFCISW.NE.1)CALL TSTERR (6,1SW) tBUS ERROR? 
D IF (CISW.NE.1) WRITE (6,110) ISW 
D116 FORMAT (' WRMEA ACKNOWLEDGE STATUS: *13) 
6010 CALL SM102 (“24,FREQ.'D °,JSW,1ISW)!DISABLE SYNTHESIZER AGAIN 
6011 RETURN 
Cc 
0012 END 


PROGRAM SECTIONS 


Name Size Attributes 
SCODE1 0060252 85 RW,1,CON,LCL 
SPDATA 900042 (7 RW.D, CON, LCL 
SIDATA 000064 26 RW,D,CON,LCL 
SVARS 060006 3 RW,D,CON, LCL 
Total Space Ailocated = 000406 131 


No FPP Instructions Generated 
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FORTRAN IV-PLUS V3.6 21:09:42 25-Mar-62 
METST .FTN;112 “TR: BLOCKS, WR 


9061 SUBROUTINE METST (IAD,1T.MODE) 


SUbhROUTINE TO TEST THE LEESBURG’ RD-11/LSI-11“A22D 
H!GU-POWER MEASUREMENT SYSTEM INPUTS. 


ACLULE 'LSPARM. COM’ 
OUMON “LSPARM/SIFREQ,STAMPL.SAFREQ.SAAMPL, 
11GAUNC?) NAD, NADT NPAD .NPTR ,NPST,NCYC,NHAR, 
2iAVG,IDST, ICTYP, IAERS, 1GATE, (WINDO,IDLUN, /HARC, 
3IPORT (I) 
SIFRE@ ,S!AM?@L----~SIGNAL FREQUENCY AND AMPLITUDE 
SAFREQ, SAAMPL SAMPLE FREQUENCY AND AMPLITUDE 

4/D PROGRAMMABLE GAIN 

WHICH A/’N’S USED (1=1, 222,274, ALL=7) 

WHICH Av} «DMA? TO INTERRUPT ON 

NUMBER OF POINTS TO CONVERT 
NPTR------NUMBER OF POINTS 10 TRANSFER 
NPST----- -FIRST POINT TO TRANSFER 

NUMBER OF YCLES TO CONVERT 

NUMBER (¢ HARMONICS "SEI TO COMPUTE DISTORTION 

--MUMBER OF CYCLES T:) VERAFE 10 CO=DON'TS 
DISPLAY STYLE (@-~%) 
---COMPUTATION TYPE 

avi ERRORS TG RETRY BEFORE QUi7T TING 
'CATE-<-~--- GATED (1) OR NCT (0) CODRV-il IN OSE)) 
UWINDO-~-~-*MULTIPLY DATA BY WINDOW (WDATA) 
\OLUN--~--DISPLAY DEVICE 
THARC HARMONIC TO COMPUTE 


AAAAASANMANAAANO 


Ot 


BRAN KS RRR RR RRR KKK HK HRRKR HRA 
~~ 


Qaaaa 


904 'NCLUDE 'LSDATA.COM’ 

HHO PARAIIETER BSIZE=1024 {MAXIMUM POINTS TO READ 

096 PAHANETER NOAD=3 {NUMBER OF CHANNELS 

BOOT COMBAN » 1 NAT/ I DATA (CNOAD*BS1ZE) 1RAW A/D DATA (IN BITS) 

00e8 COMMON ~RDAT/RDATA (NOAD*BS[ZE} 'NORMALIZED DATA (1N VOLTS) 
OBn9 CONNON /“WDAT/WDATAC(BSIZE) ,WATT .WAV !WINDOW, ATTENUATION, AVERAGE 
001¢ COUNeN “SDAT/SDATA (BSIZE) {COMPUTATION BUFFER 

Aart COMMON -ATIM. TIMCNGAD) !COMPUTATION TIMES 


RH RRR HR HK 


“Al? CNTEGER CPARCIS). TAEC6) 


VL or TS EQ. ME') GOTO 160 
ar A CIT LL RE) COTO 200 
BBLS .Y CETLEQ.* RS") GOTO 266 
WOOL C. if cty EAL"R3" + GOTO 206 


SaeMe FO Us, TG PIRFORM AP CONVERSION (ME) 

SOW TE 

WAIVE €6.108) 

muRtAT ¢' DO A/D CONVERSION’ > 
O18 iF (MODE. NE.'LO") CALL WRLVM (IAD, ‘ME’ ,@,0,1SW) 
Hoty 1h CNOSE.EQ@.°LO'» CALL WRMEA (1AD.°ME’ .SAFREQ,20@, ISW) 
VO28 IFO CUTSW.GE.@) RETURN 
Ou2t WhiVe (6,102) [Sw 
Co Pee FOSMAT oc xxx A/D CONVERSION ERROR: ‘13° *#*%') 
Hens tF (is7.EN.-100@) WRITE (6.110) 
vov4 iv ttSw. be. +101) WRITE £6,120) 
9025 FORMAT «° a-'D SEIUP FRROR t#A/D°''S OR INTERRUPT)* ) 
®O2O FORMS © ° TOO MANY POLNTS REQUESTED? > 
0027 RETURN 


READ A DATA BUFFER FROM THE LS! (A1,R2,R3) 
aelB CON TINGE 

WRITE (6,201) 

FORMAT (° READ FROM 11703") 
Dez” ANS: SECNDS(@.0) 
Poet Ls) CALL WALVM (TAD,1T,0,9, 1SW) 
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METST .. FTN; 112 “TR: BLOCKS WR 


0031 IF CISW.NE.1) WRITE (6,203) ISW 

0032 FORMAT (' WRLVM ERROR: °13) 

0033 TF (ISW.EQ.-6) J=sIBUP(1IO) 

0034 I* cTSW.NE.1) RETURN 

0035 I. «IT. EQ.°RI') J=l 

0036 IF IT. EQ@.‘'R2’) Jz2 

0037 {F .IT.EQ@.°RY’) JF3 

0038 K=(J-1) *NPTR+! 

ec9 CALL ADLSI CIAD, 1 DATA(K) ,2aNPTR,ISW) 
WRITE (6,220) (IDATA(I),1=K,.K+NPTR) 
FORMAT (° '&08) 

00406 IF (ISW.NE.1) WRITE (6,230) J,ISW 

0v41 FORMAT (° x*«x A-D DATA CHANNEL ‘[1° IN ERROR: ‘13° xxx’) 

0042 UF cISW.EQ.-6) J= I BUP(10) 

043 iF CISW.NE.1L) RETURN 


FIX TAaE DATA PROPERLY 

004-4 Eisti+¢J-fL) *NeTR 

0045 LOFLI+NPTR 

0046 DO 250 1=LI,i2 

0047 PHS rANDO°T777 IP SVACH 

0048 IF CE. GE.*4900> IM=[OR(°174000,1M) 

0e@49 ADATA‘SL)=10.08 1M ° 4008 

0056 oF CONS ENUE 

veSi TINGS) =SECNDSCRDS? 
IF (CUT. EQ.'RI').OR. (CIT. EQ.‘R2’)) WRITE (6,260) TIM(J) 
(F (CIT. £EQ@.'R3*).OR.CIT.EQ@.° RE’)? WRITE (6,260) TIM(J) 
FORMAT  «' BUFFER READ TIME='F8.3) 


0052 RETURN 
0952 END 


FROGRAM SECTIONS 
Name Size Attributes 


SCODE! 8011350 RW,1,CON, LCL 
SPDATA 6090324 RW.D,CON, LCL 
SIDATA 980050 ax RW,D,CON, LC? 
SVARS ONOuT + 36 RW,.D,.CON, LCL 
LSPARM @¢890@72 2 RW.D.OVR.GBL 
IDAT 01400% : RW,D,OVR, CBL 
ADAT 030060 RW.D,OVR,GBL 
wDAT O1OA1I8a RW,P.OVR,GBL 
~DAT O1O00G0 i RW.D,OVR, GBL 
{TIM O6001 4 RW,D,OVR. GAL 


Total Space Allocsted = @65772 193821 


APPENDIX E 


Voltmeter Program Listings 


The modules used in this program are written in PDP-1l assembly language, 
divided intc pure (ROM) code (modules 1-7) with all dynamic variables (RAM) in ! 
the final module. 


(i) Leis Command interprece:; reasurement, computation, and 
display driver 


(2) 32D A/D converter control subroutine 


(3) CONVERT Floatiag to ASCII encoding subroutine, used to prepare 
numbers for display or transfer to the host 


(4) CAPD Computation of amplitude, phase, distortion; uses DGZL 
below to obtain real and imaginary parts of voltage 


(5). DGZL Single frequency DFT subroutine: computes real and 
imaginary parts of voltage of a given spectral line by 
Goertzel's algorithm 


(6) GPTLOSR ISSE-488 bus input and output subroutines, initialize 
subroutine, and interrupt service routine 


(7) GRMS Computation by time-integration of rms, dc, positive 
neak, voltages 


(8) LSIVAR all che dynamic variables used by (1)-(6) above 


Also used are modules from DEC's RSX-1LIM System Library-~SATAN (arctangent), 
SSIN (sine), and $SQRT (square root). 


The subroutines (2)-(€7) are written tc conform with DEC standards; 
argument iists are passec through register 5 so that these routines could be 
used by other programs. Conditional assembly parameters allow dynamic 
vatiaoles to be assembi:a locally, or in an external global file such as 
LOLVAR. 


(1) MODULE LSIVM - Main program for controlling the microcomputer- 
controlled sampling voltmeter 


Internal Summary: 
(LABEL) (FUNCTION) (PAGE) 
MONITOR ROUTINES 


START Initialize system 5 


RDCODE Read string from bus; interpret type of command 5 


PARAMETER ROUTINES 


PARAM Parameter decode and load 6 
} 
FREQP Read frequency data 8 
MESDAY Read and display text string 8 
DECOD Byte string to integer decoding 9 
FPISR Floating-point processor interrupt service 9 
routine 
EXEC Interpret and direct control codes 10 
CERR Parameter 1/0 error routine ll 


CONTROL RCUTINES 


ABORT Abort process and reinitialize ll 
MESURE Make measurement 12 
WAKE Counter/display initialization 12 
PARLST Display setup parameters 13 
LFWALT Delay sub~cutine 13 
PARSND Send setup parameters to host 14 
ADSEND Send A/D buffer or window to host 14 


PRE-PROCESSING ROUTINES 


CMPUT Computation interpreter 15 
CWAVE Compute test sine wave 17 
FLOT Convert A/D buffer to floating point buffer 18 
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: oe . vases 2a nce 
a maple ee ce Sale = pyasrtintece: ermceentinlsSs 


(LABEL) (FUNCTION) (PAGE) 


WINDIN Read sequence window from host; compute 19 
average 

WINDOW Multiply data buffer by window buffer 19 

BAVER Average daca buffer to fewer cycles 20 


GDFTD Compute voltage, phase, distorticn (DFT) of 2h 

; fundamental 

: GDFT Compute voltage, ghac= (DFT) of rundanental 21 

ie 

' GSL Compute voltage phase FL) <f harmonic ie 

i TDFTD Tratsier voltage, pma- t. cistortion to host 22 
SDEFT Transfer voitzege, phase to host; display 22 
DSLD Load DFT values for display 22 
GRMSV Do statistical computaticn 23 
GPEKO Find absolute peak value sas 
~RMSV fransfer statistical data to nest; display 24 
TPER Transter peak veltage to host; display 24 
DEGAIN oorrect voltage for A/D zain 24 
DEWIND Corrteet DFT voltage ror vindow average 24 
GPWR Compute total power (second joint moment ) 25 
oJC Compute total power and ac power 26 
GMN Get rms voitage for power computation 26 
TJC Transfer total power and ac power to host 26 


DISPLAY ROUT! YES 


VMDISP Stand-alone voltmeter executive 27 
CHODLS Display, mode 9 (conversion and error count) 29 
CHIDTS Display, modes 1-3 (voltage, phase, and 30 


distortion) 
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WRITE 


PRINT 


Display, 


Display, 
phase) 


“erminal 


Terminal 


mode 4 (voltages) 


mode 5 (voltages and relative 


output (string) 


output (integer) 
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(2) MODULE A32D - Subroutine to control ADAC A/D converters 
Calling Sequence: 
CALL A32D (BUFFER, GAINS, NWORDS, ADS, GATED, CHNO, ISW) 
where the arguments are: 
BUFFER Integer array to receive contiguous sampled data 


GAINS Integer array containing gains of individual A/D 
converters 


Number of samples to digitize per channel 


A/D setup word containing which A/D's in use; which to 
interrupt from 


Integer variable specifying presence (1) or absence 
(0) of DRV11/toneburst gate 
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Integer array containing port (channel) numbers for 
each channel 


ISW Integer variable for returned status 


Internal Summary: 


(LABEL) (FUNCTION) 


SETUP ROUTINES 
A32D Argument transfer; parameter checking 


SETDMA . Load DMA bus address and word count registers; 
load pseudo-register for DMA control 


SETA2D Load pseudo-register for A/D control; 
load DMA pseudo-registers into control status 
registers 

MEASUREMENT ROUTINE 

READY Prepare for measurement 

GO Check for pulsed-system gate 
Load A/D pseudo registers into actual control 


status registers; wait for completion 
(interrupt) 


+ cana plement 


ADISR 


DMISR 


Dummy interrupt service routine 


Functional interrupt service routine 
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(3) MODULE CONVRT - Subroutine to convert a floating-point number to an ASCII 


string. 
Calling Sequence: 
CALL CONVRT (VALUE, STRING, K, L) 


where the arguments are: 


VALUE Floating-point variable to be encoded 

STRING Byte array to load encoded variable 

K Encode with K places to left of decimal point 
L Encode with L places to right of decimal point 


NOTE: Decimal point (.) and sign (+ or -) are not included in the size of K 
and L. 
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(4) MODULE DAPD ~- Subroutine to compute voltage and phase from the real and 


imaginary parts of the voltage, and to compute the harmonic 
distortion from a specific number of harmonics. 


Calling Sequence: 
CALL DAPD (X, N, K, AMP, PHASE, DSTORT, LHARM, IERR) 


where the arguments «1 


X F°--ating-point data sequence 

N Length of data sequence 

K Spectral line (number of cycles) to compute 
AMP Resulting total rms voltage 

PHASE Resulting phase in degrees 

DSTORT Harmonic distortion 

LHARM Number of harmonics used to compute distortion 
TERR Execution status 


NOTES: Ll. Lf LHARM 1, distortion is not ccmputed 

2. If DSTORT = -1, distortion was not computed 
3. If DSTORT -2, error occurred (AMP = 0) 

- .f TERR = -1, K = N (illegal condition) 


If IERR = LHARM, harmonic folds onto dec line 


it 
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(5) MODULE DGZL - Subroutine to perform a discrete Fourier transform on a 
data sequence using Goertzel's algorithm. 


Calling Sequence: 
CALL DGZL (X, N, K, XR, XI) 
wher: the arguments are: 
x Floating-point data sequence 
N Length of sequence 
Spectral line (number of cycles) to compute 
Real part of voltage 


Imaginary part of voltage 
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(6) MODULE GPIOSR - This subroutine controls input and output along the IEEE- 


GPINLI 


GPIN 


GPOUT 


488 pus. There are three primary functions: 


Initialize the interface and requisite variables; also set up 
address to use when group execute trigger (GET) is received. 


Calling Sequence: 

CALL GPINI (RESET) 
where the argument is: 

RESET Address of routine to transfer to on GET 
Listen function - read a byte string from the ILEEE-488 bus, 
terminated by a character received with EOI true, or by input 
buffer overflow. 

Calling Sequence: 

CALL GPIN (INBUF , MAXCH, NCH) 
where the arguments are: 

INBUF Byte buffer for input string 

MAXCH Maximum number of characters allowed 


NCH Number of characters received 


Talk function - set SRQ and transmit a byte string to the IEEE- 
488 bus. 


Calling Sequence: 
CALL GPOUT (OUTBUF , NCH) 
where the arguments are: 
OUTBUF Byte buffer for output string 


NCH Number of characters to output 


Also included is an interrupt service routine. The only interrupt always 
enabled is the group execute trigger (GET) interrupt, which feature is used to 
effect a program abort or restart funcricn. The bus-in and bus-out interrupts 
are enabled only when their respective subroutines are called and disabled on 
completion of transfer. 


The stand-alone driver supplied by National Instruments Corp. was not 
used for two reasons: size (keeping the total program size under 4K) and to 
implement the special abort function. However, programming this complex 
interface was not a straightforward exercise. 
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(7) MODULE GRMS - This subroutine computes parameters from statistical 
techniques. 


Calling Sequence: 
CALL GRMS (N, X, RMS, DC, PEAKMX, PEAKMN) 
where the arguments are: 
N Length of sequence 
x Floating-point data sequence 
Returned rms voltage (standard deviation) 
Returned de voltage (average) 
Peak voltage (maximum) 


Peak voltage (minimum) 


NOTE: This code was originally written by R.W. Anderson and subsequently 
modified by R.W. Luckey, both of NRL-USRD. 
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(8) LSIVAR - The following is an alphabetical list (with definitions) of the 
variables used in the software. 


for those used by the subroutines noted in parentheses. 


Exponent factor (DGZL) 

Input port for channel 1 

Input port for channel 2 

Input port for channel 3 

A/D buffers (three contiguous arrays) 
Numbers of A/D's used 

Measurement status register (A32D) 

DMA to interrupt on (ignore all others) 
A/D pseudo-register (A32D) 

Length in bytes of buffer to send to host 
Automatic retries when A/D errors occur 


Averaging mode (0 * don't, or cycles per averaged 
subsequence ) 


Buffer size allowed per channel 
12-bit integer values from A/D #1 
12-bit integer values from A/D #2 
12-bit integer values from A/D #3 
Cosine factor (DGZL) 

Cosine factor (DGZL) 

Current factor (I = CFACT*E) 
Display character buffer 

Channel number being computed 
Computation type (not used) 
Covariance (ac power) 


Start of LSIVM input parameter buffer 


All are used by LSIVM except 


DCV 
DFD 
DIST1 
DIST2 
DIST3 
DIST 
DLUN 
DSTYLE 
ECNT 
FSAM 
FSIG 
GAIN1 
GAIN2 
GAIN3 
GAIN 
GATED 
GPBIN 
GPBOUT 
GPEOF 
GPMAX 
GPNCH 


GTRIG 


Average (dc) voltage from statistical computation 
DFT data array 

Distortion of channel 1 

Distortion of channel 2 

Distortion of channel 3 

Harmonic distortion 

Display device (1 = terminal, 2 = display) 
Display style (0-5) 

Counter for A/D error retry 

Sample frequency (Hz) 

Signal frequency (Hz) 

A/D gain channel 1 (1, 2, 5, 10) 

A/D gain channel 2 (1, 2, 5, 10) 


A/D gain channel 3 (1, 2, 5, 10) 


A/D converter gains array 


Gated voltage flag (0 = no, 1 = yes) 

Input buffer address (GPIOSR) 

Output buffer address (GPIOSR) 

E-O~I found - flag (0 = no, 1 = yes) (GPIOSR) 
Maximum characters on input (GPIOSR) 
Character count for input or output (GPIOSR) 
Address to go to if “GET" received (GPIOSR) 
Number of harmonics to use to compute distortion 
Harmonics - used flag array (DAPD) 

Specific harmonic (line + 1) to compute 

Line number (DAPD) 


Bytes received from IEEE-488 bus 


Command byte 1 

Count of A/D conversions 

Count of A/D errors 

Command byte 2 (may be channel number) 
Error status of DFT 

IEEE-488 bus input packet 

Command parameter string 


LSIVM status word 


What DMA to interrupt from (A32D) 


Pointer to displayed parameters 

Counter for displayed parameters 

Parameter to display 

Index for displayed parameters 

Sequences to average (NCYC/AVGCOD) 

Cycles per sequence 

Cycles to compute (AVGCOD); averaging changes 
Cycles to compute (for a harmonic) 

Samples per sequence (DAPD) 

Samples to convert per DMA transfer per channel 
Phase from DFT 

Peak (-) voltage from inspection 

Peak (+) voltage from inspection 

Phase of channel 1 

Phase of channel 2 

Phase of channel 3 

Total power (ac and dc) 


Phase sign (DAPD) 


PTC 


STACK 
STPT 
STRING 
Tl 

T2 


T3 


Points to compute (PTR/NCYC); averaging changes this 
Samples to transfer to host (samples per sequence) 
Harmonic power (DAPD) 

Argument list for writing buffer to host 

Start of program RAM space 

Address of buffer to send to host 

One voltage buffer (converted to floating-point) 
Statistical data array 

RMS voltage in channel 1 

RMS voltage in channel 2 

RMS voltage in channel 3 

RMS voltage from statistical analysis 

Phase difference between channels 1 and 2 

Sine factor (DGZL) 

Second joint moment (total power) 

Start (éop) of program stack 

Number of first sample to transfer to host 

Used by print 

DMA #1 pseudo-register (A32D) 

DMA #2 pseudo-register (A32D) 

DMA #3 pseudo-register (A32D) 

Total rms voltage from DFT 

What A/D converters used (A32D) 

Average wincow magnitude 

Window buffer (floating-point) 

Window mode (0 = no, 1 = yes) 


Imaginary voltage part (DAPD) 


XREAL Real voltage part (DAPu) 
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APPENDIX F 


Voltmeter Memory Map 


oc0000 
| ROM | PROGRAM CODE (4 KW) 


020000 


PRESENT PROGRAM VARIABLE SPACE (8 KW) 


060000 
‘ MEMORY AVAILABLE FOR EXPANSION (8 KW) 
120000 
ADDRESS SPACE AVAILABLE FOR EXPANSION (12 KW) 
177777 


Fig. Fl - Voltmeter memory map. 


The ROM code used by this version of the voltmeter occupies approximately 
15300 (octal) bytes of memory, so there is room for expansion. 


The program variable space allows, at present, data buffers of 1024 
samples, but this can be expanded by modifying the parameter BSIZE in LSIVAR 
when the program is rebuilt. Note that for each sample a buffer is expanded, 
memory is used as follows: 


6 bytes for the sampled A/D waveform (2 bytes per channel), 
4 bytes for the integer to floating-point conversion biffer, 
4 bytes for the floating-point data window. 


This requires approximately 7K words of RAM memory for each additionel 
1024 samples. 
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